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ABSTRACT 

Dark- matter haloes are supposed to form at the positions of maxima in the initial matter density field. The gravitational-shear field's 
ellipticity and prolaticity that serve as input for the ellipsoidal-collapse model, however, are derived from a distribution that does not 
take the additional maximum constraint into account. In this article, I quantify the variations of the most probable and the expected 
values of the ellipticity and the prolaticity when considering this additional constraint. Based on the statistics of Gaussian random 
fields, it is possible to set up a joint distribution for the eigenvalues of the gravitational-shear tensor and the matter density that 
incorporates the maximum constraint by invoking a vanishing first derivative and a negative definite second derivative of the density 
field into the calculation. Both the most probable and the expected value of the ellipticity calculated under the assumption of a 
density maximum are about 1 - 8 % lower compared to the standard distribution used in the literature in the density range relevant for 
cosmological structure formation. Additionally, the analogous quantities for the prolaticity do not vanish but acquire slight positive 
values in the range of 1CT 7 to 1CT 2 . For large overdensities, the predictions from both distributions converge. The non-vanishing 
prolaticity may have consequences for the ellipsoidal-collapse model especially when the linear model for the external shear is used. 
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O 1. Introduction 

<Jl Galaxies and galaxy clusters as collapsed objects are supposed to 
, ^ , form from peaks in the initial matter density field. After an over- 
dense patch has decoupled from the Hubble expansion, the ex- 
i—l tent of this patch is supposed to reach a maximum value, then to 
J> shrink again until the collapse stops when the structure reaches 
I"*"- virial equilibrium. Analytical models for this evolution include 
the spherical-collapse model (e.g. Wang & Steinhardt 1998; Pace 
"^t* et al. 2010) and the ellipsoidal-collapse model (e.g. Bartelmann 
O et al. 1993; Eisenstein & Loeb 1995; Bond & Myers 1996). 

From the theory of inflation, the density field is supposed to 
< O l be a random field whose non-Gaussianity is small (Planck Col- 
laboration XXIV 2013). Since also derivatives and integrals of 
\ I Gaussian random fields are again Gaussian, a multivariate Gaus- 
J>. sian random field can be used to derive the initial number den- 
: T ~j sity of peaks in the density field as a function of peak height 
by applying respective constraints to the density field's first and 
second derivatives (Bardeen et al. 1986). 

Not only can the density field be described as a Gaussian 
random field, but through Poisson's equation, also the gravita- 
tional field has to follow a Gaussian distribution and therefore 
also the eigenvalues of the gravitational-shear tensor can be de- 
rived from the statistics of Gaussian random fields (Zel'dovich 
1970; Doroshkevich 1970). These eigenvalues can be combined 
to define the ellipticity and the prolaticity of the gravitational- 
shear field whose distribution for a given value of the overden- 
sity can be derived by transformation of variables (Sheth et al. 
2001) and either their most probable or their expected values are 
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important ingredients for the ellipsoidal-collapse model by Bond 
& Myers (1996). 

However, this distribution does not take into account that 
haloes should form at maxima in the density field, but rather 
yields the probability density for ellipticity and prolaticity for 
any overdensity irrespective if it is taken at a maximum or not. 
Recently, attempts have been made to take the maximum con- 
straint into account when calculating the gravitational-shear ten- 
sor's constrained eigenvalues (Rossi 2012) as well as the con- 
strained ellipticity and prolaticity from them (Rossi 2013) by ex- 
tending the formalism of Doroshkevich (1970) and Sheth et al. 
(2001) in the way to also include the density's Hessian. One cru- 
cial condition for a maximum, namely the vanishing first deriva- 
tive of the density, however, is not taken into account. Addition- 
ally, it does not seem possible in the framework of Rossi (2012, 
2013) to derive the distribution for the unconstrained ellipticity 
and prolaticity necessary for the ellipsoidal-collapse model since 
the gravitational-shear tensor and Hessian of the density cannot 
be diagonalised simultaneously. 

In this article, I will show how the distribution of the gravita- 
tional shear's ellipticity and prolaticity at the position of density 
maxima can be derived from the statistics of Gaussian random 
fields by accounting for both the vanishing first derivative of the 
density and the constraint that the density's Hessian has to be 
negative definite. In § 2, 1 give a short introduction to Gaussian 
random fields and show step-by-step how the targeted distribu- 
tion can be calculated invoking the technique to derive number 
densities from multivariate Gaussians of the relevant quantities 
based on Bardeen et al. (1986) and how the constraint for the 
density's Hessian can be properly considered without diagonali- 
sation. In § 3, 1 compare the results from the distribution includ- 
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ing the maximum constraint to the one by Sheth et al. (2001) 
that does not take it into account. Finally, I summarise the main 
results of this article in § 4 and briefly hint at implications for 
the ellipsoidal-collapse model by Bond & Myers (1996). 

2. The conditional distribution for ellipticity and 
prolaticity 

In this section, I derive the conditional probability for the ellip- 
ticity and the prolaticity of the gravitational shear from the statis- 
tics of Gaussian random fields. First, I give a brief overview over 
the theoretical background of Gaussian random fields based on 
Bardeen et al. (1986) (see also Angrick & Bartelmann 2009 for 
further information). Hence I also adopt F(r) for the random 
field and T)(r) = Vf(r) and £,y(r) = didjF{r) for its first and 
second derivatives, respectively. 

2. 1 . Theoretical background 

An n-dimensional random field F(r) assigns a set of random 
numbers to each point in n-dimensional space. A joint proba- 
bility function can be declared for m arbitrary points r ; as the 
probability that the field F, considered at the points r ; , has val- 
ues between F(jrj) and F(rj) + dF(ry) with 1 < j < m. 

A Gaussian random field is a field whose joint probability 
functions are multivariate Gaussians. Let with 1 < i < p 
be a set of Gaussian random variables with means {yi) and 
Ai/j = yi — (yi). The covariance matrix M has the elements 
Mjj = (At/,- Ayj), and the joint probability function of the Gaus- 
sian random variables is 

1 

: e ^dy x - --dyp (1) 



p(y\ 



> y P ) dyi ■ ■ ■ dy P 



V(27t)''detM 



with the quadratic form 
g^I^A^M- 1 ) Ayj. 



(2) 



In the following, I will only consider fields on the 3-dimensional 
Euclidean space with coordinates r = {x\,X2,x^) l: . A homo- 
geneous Gaussian random field with zero mean is fully charac- 
terised by its two-point correlation function ^ F (r 1 ,r2) = ^p(\ r \ - 
T2\) = (F(ri) F(r2)) or equivalently its Fourier transform, the 
power spectrum Pf(fe). Its spectral moments are defined by 

dk P F (k)k^ 2 W 2 R (k), (3) 



°7 = 



2n 2 



where Wr(Ic) is the Fourier transform of the filter function W«(r), 
and R is its filter size. In the following, I adopt Gaussian filter- 
ing, hence V^(r) = exp[-|r| 2 /(2/? 2 )] and W(k) = exp(-R 2 k 2 /2). 
The zeroth-order spectral moment cr 2 is the power spectrum's 
variance. 

It is possible to scale the field F as well as its derivatives 
with corresponding spectral moments so that they become di- 
mensionless. I will denote these reduced fields with a tilde, e.g. 
F = F/cTo, fj = t\Ict\, and £ = C/o"2- A quantity that will be 
used more often in the next sections is a special combination of 
the first three spectral moments, y = <x 2 /(cr o~i). 

The elements ip,j of the gravitational-shear tensor cp 
(Zel'dovich 1970) contain the scaled second derivatives of the 
gravitational potential 

1 

yi, - (4) 



where G is Newton's constant and pb the background matter den- 
sity, so that tr (p = 8 by Poisson's equation. Here, 6 is the dimen- 
sionless overdensity defined by S = {p-pb)lpb- I will denote the 
eigenvalues of cp with A/, where 1 < i < 3. The ellipticity e and 
the prolaticity p are combinations of these eigenvalues, 
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M - 2A 2 + A3 
25 



with 5 - A1+A2+A3. (5) 



Ordering the eigenvalues in the way A\ > A2 > A3, as it is com- 
monly done in the literature, imposes the constraints e > and 
—e < p < e for 6 > 0. Since a spherical configuration has 
Ai = A2 = A3, this is equivalent to e = p — 0. While the ellip- 
ticity measures the anisotropy in the xi-X3-plane, the prolaticity 
quantifies the shape perpendicular to it: if p < 0, the configura- 
tion is prolate, for p > it is oblate. Note that ellipticity and 
prolaticity are defined on the linear density field and should not 
be confused with the shape of the final halo! 

Following Bardeen et al. (1986), the conditional probability 
for the random variables e and p given a maximum of reduced 
height 6 can be expressed as 



P(e, p 1 6, max) 



J ad d 6 l\detl\P(e,p,S,f\ =0,1) 
f Bi d 6 l\detZ\P(5,j\=0,Q 



(6) 



where here and in the following, fj is the reduced first deriva- 
tive of 5 and C is the corresponding Hessian, which has only six 
independent components since it is symmetric. Note that only 
those subvolumes contribute to the integral for which £ is neg- 
ative definite (indicated by n.d.), i.e. its three eigenvalues are 
negative. How this can be properly taken into account, I will 
discuss in Sect. 2.3. 

2.2. The joint probability of relevant quantities 

To calculate P(e, p \ 6, max), I shall start with establishing 
the 15 -dimensional joint probability of reduced variables 
f(<p, fj, C) d 6 if d 3 ?7 d 6 £ to evaluate the numerator of Eq. (6). 

The correlation matrix M incorporates all auto- and cross- 
correlations of the various random variables. These are 



(ipijipki) = y^r Kiju, {flifjj) = i 8tj, 

„ „ 1 y 

{£,ij£,kl> = Kijkl, (<Pij£kl) = - j-g K ijkl< 



(7) 



where <5y is the Kronecker symbol and the quantity 
Kjjki is a combination of Kronecker symbols given by 
Kiju = djjdki + 5ikdji + Sudjk- All other correlations apart 
from those listed above vanish (cf. also Bardeen et al. 1986; 
Rossi 2012). Defining the vector of random variables y = 

(011,022,^33, 012,013,023, ^1,^2,^3, ?11,?22,?33, ?12, ?13, ?23) T , 

the correlation matrix M looks like 



M = 



1 
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where the different submatrices are given by 
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The determinant of M is 

detM= V. V , 
3 i5. 5 io 

and the inverse of M is 
M" 1 = 



2(1 - r 2 ) 
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2.3. Sylvester's criterion 

To evaluate the numerator of Eq. (6), I have to integrate over 

(10) those parts of the volume spanned by the matrix elements of C 
for which the latter is negative definite. This can be achieved by 
taking into account the following criterion. A Hermitian matrix 
is positive definite if and only if all its leading principal minors 
are positive, where the i-th leading principal minor is the de- 
terminant of its upper left i x i sub-matrix. This is known as 
Sylvester's criterion. 

Since the negative of a positive definite and Hermitian ma- 
trix is negative definite, this implies that all odd principal minors 
have to be negative and all even principal minors have to be posi- 
tive. Technically, when integrating Eq. (15) numerically over the 

(11) various matrix elements, I multiply it by a function /(£,) defined 
as 



where 1 and are defined as above and 
C 



( 4 -1 -\\ 
1 4 -1 
1 -1 4 



(12) /(O = 



Rotating the coordinate system such that cp becomes diagonal 
and assuming that the eigenvalues are ordered like ~A\ > A2 > A3, 
Bardeen et al. (1986) showed that the volume element d 6 ijo can 
be written as 



d 6 ip = 2tt 2 (2 1 - ~A 2 ){Ay - A 3 )(A 2 - A 3 ) d\ d~A 2 dA 3 . 



(13) 



Note that the two matrices qj and £ do not commute and hence 
cannot be diagonalised simultaneously. Therefore, I still have to 
keep the off-diagonal elements of t in the following calculation, 
whereas <p\2 — ipn = $23 = and if a - A, after the diagonalisa- 
tion. 

Switching from the eigenvalues to the new variables 6, e, 
and p, the Jacobi determinant of the transformation introduces 
an additional factor so that 



1 iffn<0A 



else, 



£ll ?12 
<Tl2 ?22 
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In 
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£33 



<o, 



(18) 



and integrate over all six elements from -00 to 00. 



2.4. The number density of maxima in the density field 

Analogous to the way I have derived Eq. (15), it is possible to 
derive the joint probability P(6, fj, C)d5d 3 ^d 6 |'. A similar and 
more detailed description of the various steps can be found in 
Angrick & Bartelmann (2009), where I have derived the number 
density of minima in the gravitational potential. Rotating the co- 
ordinate system such that t becomes diagonal and transforming 
to the variables 6^, e^, and p% introduced in Eq. (16), which are 
combinations of the eigenvalues X^j of C after the rotation, the 
integrand in the denominator of Eq. (6) can be written as 



Z ~9 ~ 

dA\dA 2 dA 3 = -6 dddedp. 



(14) 



Altogether, the integrand in the numerator of Eq. (6) can be writ- 
ten as 



P(e, p, 6, r\ = 0, 1 det £| 



P(S,r\ =0,Q|detCI 
25 Vl5 



57T 3 yl - j' 



x exp 



e ( {e} - p})\[(\ + P( ) 2 - 9e 2 ](2 P( - 1)| 8* 
(3 2 + 3 2 + 2y55 ( + (1 - y 2 )(15e 2 + 5 P( )S 2 



2(y 2 - 1) 



(19) 



.13/2 



e ( e 2 _p 2 ) g 5 eX p J - [3(2<5 2 - 5A) where I have ordered the eigenvalues of C like A^\ 
\2{y 2 - 1) To evaluate the denominator of Eq. (6), I have to i 



,2) n n > 2 (l -y 2 ) 3 
+ 2ySS ( (l + I5ee c + 5pp ( ) + 6 2 (l + 15e 2 + 5p 2 )] \ |det l\ 



(15) 



> A^ 2 ^ ^,3- 

To evaluate the denominator of Eq. (6), I have to integrate over 
the eigenvalues A ( j from -00 to 0, or, since I have transformed 
to the new variables S^, e^, and pr, integrate over them in the 
proper range. The condition > A^\ > A^ 2 ^ \,3 translated to 
the new variables reads 



where I have introduced the abbreviations 



A = £n £22 + luh-i + lnhi - In ~ In ~ Hi' 



6 ( < A 



c _ y , Z , ~y _ In ~ha _lu -2f 2 2 + ?33 

°( = £11 + in + £33, e ( = — ^= , p ( = - 



25r 



26r 



~ < e c < ~ A -1 - 3e ( < p ( < -e ( 



V I — — < e% < A e% < p( < —en 



(20) 



(16) 



and with 



det I = lutein + llulnU ~ Inln ~ luln ~ hslli- (17) 



Note that the variables 6^ and e^ are negative in contrast to the 
analogous variables 3 and e for the gravitational-shear tensor. 

It turns out that the integrations over e^ and p^ can be done 
analytically, while the integration over has to be done numeri- 
cally. The denominator of Eq. (6) is equal to the number density 
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v = 1 



v = 5 




of maxima n maK (6) of the reduced density field 6 if multiplied by 
a factor cr^/cr\. The latter can be written as 



«max(5) = 

with 
// = exp 



of ro 

o\ J-oo 



A5r 



12007t 3 V3(l - y 2 ) 



(V5y + 25 V27tw) 



(21) 



6 2 +2y66 ( + S 2 (6-5y 2 ) 



2(r 2 - 1) 



v = 106j - 32 + (155<J 2 + 32) exp 



(15$ 2 \ 



w = 6( (6 2 - 3) exp 



erf 



erf 



(22) 



Now I have all the ingredients to express the probability to 
find the ellipticity e and the prolaticity p of the gravitational- 
shear field in the range [e, e + de] and [p, p + dp] given a density 
peak of height 6. To summarise, after scaling it to the reduced 
density 6 = 6/<To, it is given by the ratio (6), where the numer- 
ator is expressed by Eq. (15) together with Eq. (16), and the de- 
nominator by Eqs. (21) and (22) except the factor o-\lo-\. Equa- 
tion (18) offers a nice way to take into account the constraint 
that Eq. (15) has to be integrated only over those subvolumes for 
which the matrix C is negative definite. 



2.5. Comparison to the ansatz of Rossi 

Rossi (2012) derives a conditional probability for the reduced 
shear tensor cp given positive definite Hessian C (see his Eqs. (18) 
and (19)). In the language of this paper, his probability looks like 



and the matrix C to use their eigenvalues in the further calcu- 
lation. Additionally, both P((p \ Q and P{C) can be written as 
an expression analogous to the distribution of the gravitational- 
shear tensor's eigenvalues derived by Doroshkevich (1970). 

However, it is not possible to diagonalise both matrices <p | C 
and C simultaneously since the commutator of the two does not 
vanish. Hence, the aim of Rossi (2013) that this ansatz can be 
used to determine the probability P(e, p \ 6, max) from Eq. (23), 
implying that the eigenvalues of both the former matrices can be 
used, is not possible. Since only one of them can be diagonalised 
at once, it is not obvious how the six-dimensional integration 
over all entries of the matrix £ can be circumvented. 

Furthermore, Eq. (23) is incomplete in the sense that it ne- 
glects the constraint that the first derivatives r| have to vanish 
for a maximum. Rossi (2012, 2013) only take into account the 
Hessian which has to be positive definite (actually, it should be 
negative definite for a maximum) and therefore their probabil- 
ity is not only based on maxima in the density field but also on 
areas for which the first derivative is non-zero. Equation (6) is 
more specific in that respect since it takes the constraint r| = 
explicitly into account (see also Bardeen et al. 1986). 



3. Results 

In this section, I present some numerical results based on the 
theoretical work developed in the previous section. For the 
calculation, I adopt a flat reference ACDM model with mat- 
ter density £2 m = 0.3, dark-energy density Q. A = 0.7, baryon 
density Ob - 0.04, a Hubble constant of h — 0.7 in units of 
100 km s _1 Mpc~', and power spectrum normalisation <x 8 = 0.8. 
All previous values are set at redshift z — 0. The primordial 
spectral index of the power spectrum is set to n s — 1, and I adopt 
the transfer function of Bardeen et al. (1986). 

In the following, I compare the distribution (6) with the dis- 
tribution 



P(q>|Cp.d.) = 



f pd d 6 ZPCQP«p\~Q 



J P .d. d6 W) 



with P(<p 1 1) = 



P(<P, P(e,p\6) 
P(l) ' 



1125 



p 2 )S 5 exp 



5« 



-6\3e z + p L ) 



(24) 



(23) 



where p.d. indicates integration over the subvolume for which C 
is positive definite. The reason for this complicated way of writ- 
ing is the aim to diagonalise both the conditional matrix <p | C 



(Sheth et al. 2001), which expresses the same probability as 
Eq. (6) except that the maximum constraint is not taken into ac- 
count, for different values of 



y = 6 2 Jal(R). 



(25) 
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Fig. 2. The distributions P(e, p | v, max) and P(e, p | v) for fixed values of either e or p. Left panel: p = for v = 1 and v = 5. Central panel: 
e = 0.2 for v = 1 and v = 5. Right panel: e = 0.7 for v = 25. 



Here, 6 C is the linear overdensity of the spherical-collapse model, 
and o-()(R) is the power spectrum's variance exponentially fil- 
tered on the scale R (cf. Eq. 3). Hence, for a given v, the radius 
R has to be chosen such that Eq. (25) is fulfilled. 

In Fig. 1, 1 compare the distributions (6) and (24) as a func- 
tion of e and p for four different values of v. For each v, the 
corresponding reduced density contrast is simply 6 = -\/v. From 
Eq. (25), one can infer a corresponding filter radius R that is then 
used to calculate the higher spectral moments (Ti(R) and <T2(R) 
from Eq. (3) and then also y(R) = <t\{R)I[o-q(R) cr 2 (fl)]. 

The difference between the two distributions P(e, p | v, max) 
and P(e,p | v) is the larger the smaller v is. Taking the maxi- 
mum constraint into account, the distribution for e and p peaks 
at slightly lower values of e, implying that haloes of lower mass 
tend to be slightly less elliptical than predicted by the standard 
formula (24). As a consequence of the normalisation to unity, 
also the peak's height is larger. However, the higher v is - and 
therefore the larger a halo's mass is - the more the difference 
between the two distributions decreases. 

The physical reason why the two distributions should be 
more similar for large values of v instead of small ones is the 
following. For large v, the Gaussian random field is filtered over 
a large volume, and maxima at 6 — 6 C correspond to maxima 
whose heights are times the standard deviation. The num- 
ber density of points with v = 25 corresponding to 5 c Ictq - 5 
is fairly small and the probability that such a point is addition- 
ally a maximum is fairly high. Hence, the requirement that these 
points correspond to maxima in the density field is almost al- 
ways fulfilled naturally and then does not impose an additional 
constraint. In the limit v — > oo, the two distributions should be 
identical. 

In Fig. 2, I present three cuts through the distributions 
P(e, p | v, max) and P(e,p\v) to highlight the differences be- 
tween them. The cut at vanishing prolaticity shown in the left 
panel for v = 1 and v = 5 underlines that the former peaks at 
lower ellipticity compared to the latter. 

The central panel presents a cut at e — 0.2 again for v = 1 and 
v — 5. Contrarily to P(e,p\v), the distribution P(e, p | v, max) 
does not peak at p = for a given e but is shifted towards small 
positive values for p, which can be only seen by eye for v = 5. 

To stress this additional feature which is not obvious from the 
2D distributions in Fig. 1, 1 present a cut for e = 0.7 at v = 25 in 
the right panel of Fig. 2. At this high value for v, the amplitudes 
of both distributions are very low (~ 10~ 80 for P(e, p | v, max) 
and ~ 10~ 70 for P(e,p | v)) but the asymmetry of the former is 
now clearly visible: it peaks at p « 0.13 instead of 0. 

What is most important for practical purposes, however, is 
not the full distribution for e and p, but only the values that are 
either expected or most probable. The expectation values of e 



and p are used e.g. by Angrick & Bartelmann (2009) to deter- 
mine the initial eigenvalues for the ellipsoidal-collapse model 
by Bond & Myers (1996), whereas the most probable values are 
used by Sheth et al. (2001) to determine the ellipsoidal barrier 
when establishing their mass function. The expectation values 
for e and p as well as their variances can be calculated analyti- 
cally for the distribution P(e,p \ v). They are as a function of v 

3 2 19tt-54 

VIOttv 60ttv 

<P> = 0. 0* = ^ (26) 

(Angrick & Bartelmann 2009). The most probable values for e 
and p from the same distribution are 

V5v 

respectively (Sheth et al. 2001). 

In Fig. 3, I present the expectation values for both elliptic- 
ity and prolaticity as well as their 1-cr regions for the two dis- 
tributions P(e, p | v, max) and P(e, p\v) as a function of v. As 
expected, the expectation value for the ellipticity is smaller if 
the maximum constraint is taken into account and the expected 
prolaticity has a small non-vanishing value. For larger v, the dif- 
ference between both distributions decreases. While the largest 
difference in (e) is at v m 1 .4 with a value that is ~ 8 % lower 
if the minimum constraint is taken into account, the difference 
in (p) is monotonically decreasing with increasing v and has a 
value of ~ 0.0085 at v = 0.1 and ~ 0.005 at v = 1.4. For v = 25, 
the difference in (e) has dropped to ~ 1.3 % and ~ 0.00015 for 
(p). At the borders of the 1-cr regions, the discrepancy between 
the results from the two distributions (6) and (24) is of equal 
magnitude as for the expectation values themselves. 

Although P(e, p \ v, max) peaks at relatively large values of p 
for large e and v as shown in the right panel of Fig. 2, its ampli- 
tude has already fallen off by several orders of magnitude com- 
pared to the maximum so that this effect does not have any large 
influence when calculating the expectation value of p. Hence, 
(p) calculated from (6) does in fact approach for v — > oo rela- 
tively quickly. 

In Fig. 4, I compare the most probable values e mp and p mp 
from the two distributions with each other. Qualitatively, the 
dependence on v is like for the most probable values for both 
ellipticity and prolaticity. At v w 1.5 the relative difference in 
e mp is maximal with a value that is ~ 8 % lower if the minimum 
constraint is taken into account. At v = 25 the difference has 
dropped to 1.4 %. The most probable value for p is ~ 0.006 for 
V = 0.1 and has dropped to ~ 0.00005 for v = 25. 
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Fig. 3. Expectation values for the ellipticity (left panel) and the prolaticity (right panel) as well as their 1-cr regions for the two distributions 
P(e, p | v, max) and P(e, p | v) as a function of v. Quantities derived from the two distributions are denoted by the subscript "max" and "Sheth", 
respectively. At the bottom of each plot, the differences between the results from both distributions are quantified. In the left panel, the relative 
difference is given by A(e) lc i = ((e)sheth - (e>max)/(e>sheth ! whereas in the right panel, the absolute difference is simply A<p) a b s = </?)max - (p)sheth = 

\P>max- 



In summary, the ellipticity used in analytical models of struc- 
ture formation is about 1 - 8 % too low in the range of v relevant 
for cosmology since the maximum constraint is not taken into 
account. The prolaticity is usually assumed to vanish, but should 
be slightly positive with values at the level of 10~ 5 to 10~ 2 in the 
same range. 



4. Summary & conclusions 

In this article, I have derived the conditional probability for the 
gravitational-shear field's ellipticity and prolaticity at the posi- 
tion of a peak in the density field with a given height. The deriva- 
tion is based on the statistics of Gaussian random fields (Bardeen 
et al. 1986) and extends the works by Sheth et al. (2001) and 
Rossi (2012, 2013) as follows. In contrast to Sheth et al. (2001), 
I have added explicitly the maximum constraint in the calcula- 
tion since it is expected that haloes form at the positions of peaks 
in the density field. Although Rossi (2012, 2013) started go- 
ing into the same direction, they have considered the conditional 
eigenvalues of the matrix <p 1 1 and the corresponding elliptic- 
ity and prolaticity constructed from the constrained eigenvalues 
but not the unconditional ones that serve as ingredients for the 
ellipsoidal-collapse model by Bond & Myers (1996). Addition- 
ally, they only considered that the density's Hessian has to be 
negative definite for a maximum but did not take into account 
that the first derivative r| has also to vanish. 

The main results of this article can be summarised as follows: 

- The gravitational-shear tensor cp and the density's Hessian 
cannot be diagonalised simultaneously. Since the eigenval- 
ues of cp are needed for the definition of the ellipticity e and 
prolaticity p, the priority is on the diagonalisation of cp so 
that first, the constraint that the density's Hessian C has to 
be negative definite for a maximum can only be considered 
via Sylvester's criterion and second, the full six-dimensional 
integration over d 6 £ has to be carried out. 

- The conditional probability for e and p given a density max- 
imum of height 6 is given by (6) after rescaling to 6 = 6/o-q. 
The numerator can be calculated by Eqs. (15) and (16), 
while the denominator is given by Eqs. (21) and (22). Equa- 
tion (18) takes into account that one has to integrate Eq. (15) 
only over those subvolumes for which the Hessian C is neg- 
ative definite. 



- The distribution of e and p incorporating the maximum con- 
straint peaks at lower ellipticities and slightly positive pro- 
laticities compared to the unconstrained one by Sheth et al. 
(2001). 

- As a consequence, both the expected and the most proba- 
ble ellipticity are about 1 - 8 % smaller for 0.1 < v < 25 
with the largest deviation at around v ~ 1.4-1.5. Both the 
expected and the most probable prolaticity have slightly pos- 
itive values at the level of 10~ 5 to 10~ 2 with increasing value 
for decreasing v. 

- Although the most probable prolaticity can have relatively 
large values for large v and e, the probability density is sev- 
eral orders of magnitude lower compared to the peak value 
so that this yields only moderate changes in {p) and p mp . 

- Both expected and most probable values for ellipticity and 
prolaticity from the distribution P(e,p\ v, max) tend to con- 
verge to the values from the distribution P(e, p | v) for v — > oo 
since the larger v, the more often is the maximum constraint 
automatically fulfilled. 

The additional maximum constraint alters the main results 
from the distribution for e and p only moderately and especially 
for small values of v. Massive haloes should therefore be less 
affected from this conclusion. 

However, it will be interesting to see if the slightly positive 
values for (p) and p mp do have an impact on the ellipsoidal- 
collapse model by Bond & Myers (1996) when using the lin- 
ear approximation for the external shear since in this model, the 
eigenvalues of the external shear tensor A ext j scale with the lin- 
ear growth factor of structure formation with time. From the 
P(e, p\v), it follows that /l ex t,3 = _ ^ext,i and /t ex t,2 = so that 
^ext,2 = stays zero over the whole period of time when the 
ellipsoid first expands and then collapses. This, however, will 
change if /t ex t,2 + for small times since it will grow over time 
according to the linear growth factor of structure formation. 
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Fig. 4. Most probable values for the ellipticity {left panel) and the prolaticity (right panel) for the two distributions P(e, p \ v, max) and P(e, p \ v) 
as a function of v. Quantities derived from the two distributions are denoted by the subscript "max" and "Sheth", respectively. At the bottom of 
the left plot, the relative difference between the results from both distributions is quantified. Further denotations are as in Fig. 3. (The small bumps 
in the curves, especially in the right panel, occur from numerical noise.) 
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